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Lepton asymmetries from neutrino oscillations 

R. R. Volkas^* 

'^School of Physics, Research Centre for High Energy Physics, 
The University of Melbourne, Victoria 3010, Australia 

Reasonably large relic neutrino asymmetries can be generated by active-sterile neutrino oscillations. After 
briefly discussing possible applications, I describe the Quantum Kinetic Equation formalism used to compute 
the asymmetry growth curves. I then show how the basic features of these curves can be understood on the 
basis of the adiabatic limit approximation in the collision dominated epoch, and the pure MSW effect at lower 
temperatures. 
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1. BRIEF OVERVIEW 

The a flavour relic neutrino asymmetry L^^ 
defined by 



IS 



(1) 



where n^ is the number density of species ip, and 
a = e, /J,, T. Reasonably large values for L^^ will 
be generated by Va <-^ i^s and i7a ^^ Vs oscillations 
(where s denotes a sterile species), in the temper- 
ature range 1 MeV to lO's of MeV, provided that 
Am^ < and the vacuum mixing angle 9q is small 
jl]-D. (Note that |Ato^| must be sufficiently large 
to induce significant asymmetry growth prior to 
neutrino decoupling at 1 MeV. We focus on this 
case.) By convention, Arn? < means that the 
predominantly sterile mass eigenstate is lighter 
than the predominantly active eigenstate. The 
asymmetry growth is seeded by the CP asymmet- 
ric background plasma (baryon/electron asymme- 
try), and it is driven by runaway positive feedback 
conditions that prevail for a short time below a 
critical temperature Tc. Figure^ shows \L,^^\ evo- 
lution for three different vacuum oscillation pa- 
rameter choices M. 

It is noteworthy that this mechanism requires 
light sterile neutrinos. These hypothetical de- 
grees of freedom are interesting for other reasons 
also: as a simple explanation for two-fold maxi- 
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mal mixing (via the mirror matter [g| or pseudo- 
Dirac schemes Q), and in order to reconcile the 
solar 01 , atmospheric |g] and LSND anomalies . 
Large neutrino asymmetries have found three 
potential applications. The first is the suppres- 
sion of sterile neutrino production |l|,p|,p|JlC|| . Re- 
call that the matter-affected mixing angle for neu- 
trinos is given by 



sin^ 29„,, = 



(b 



(2) 



where s = sin20o and c = cos20o- The functions 
a and b contribute to the effective potential [O 



Va 



At 



2p 



-{-a + < 



(3) 



where p is the neutrino momentum or energy and 
_ 2pV2GFn^L 



b = 



2x/2p2 Gf n^ 



AaT 



Am 



(4) 



•-w 



The quantity Gp is the Fermi constant, m^i/ is 



the W boson mass, A^, ~ 17, A 
effective total lepton number L is 



^,T 



4.9 and the 



L = L^^+ L^^ + L^^ + L^^ + T], (5) 

where |?7| ^ 5 x 10^^° is a small term due to the 
asymmetry of the electrons and nucleons. The a 
term is the Wolfenstein effective potential, while b 
is a finite temperature correction that is very im- 
portant in the early universe context. When Va 
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Figure 1. Examples of lepton asymmetry growth curves driven by i^r <-> J^s and the corresponding 
antineutrino oscillations. The mixing angle is selected to be sin^ 26 = 10^*. The three curves correspond 
to Am^ = —0.5, —50, —5000 eV^, reading left to right. This figure is taken from Ref.B. 



is large, the matter-affected mixing angle is small 
and sterile neutrino production is suppressed. At 
high T, the b term alone is sufhcient for this pur- 
pose, but at low T a large enough L is required. 
The second application requires an electron- 
like asymmetry L^^ of the appropriate magni- 
tude. Primordial light element abundances de- 
pend on the expansion rate of the universe dur- 
ing Big Bang Nucleosynthesis (BBN) as well as 
the reaction rates for V^p <-^ e+n and Ve^ ■*-*■ e~p. 
The latter are affected by L^^ , with a positive 
asymmetry lowering n/p and thus also the He- 
lium abundance. Now, the Helium yield in- 
creases with the baryon density tjb- A posi- 
tive L^^ thus allows one to obtain an acceptable 
"^He mass fraction with a higher baryon density 
than would otherwise be the case |g|,0. Inter- 
estingly, the recent Boomerang/Maxima cosmic 
microwave background anisotropy measurements 
prefer a relatively high rjB Jl^]. This is a hint 
that Li,^ > should be taken seriously for phe- 
nomenological reasons. 



The final application is related to the second. 
Recall that the effective potential depends on the 
linear combination L which contains the baryon 
asymmetry r/. Because of this, neutrino asym- 
metry generation depends on the properties of 
the baryonic component of the plasma. If the 
latter were to be spatially inhomogeneous, then 
it could seed inhomogeneous asymmetry creation 
[nj]. There could be a domain structure, whereby 
the sign of L^,^ flips between spatial cells. This 
will become especially interesting if future obser- 
vations uncover firm evidence for inhomogeneous 
light element abundances. 

Another interesting issue is the possibility of 
rapid oscillations in the asymmetry immediately 
after Tc [|l|, [l5|jlq ] . Recent work has shown that 
there are definitely no L^^ oscillations for a large 
region of oscillation parameter space ||l6| . The re- 
maining region requires significant computational 
effort, and has not been properly explored. There 
are indications for rapid oscillations, but numeri- 
cal error cannot be definitively ruled out as their 



cause. On the theoretical front, it is strongly sus- 
pected that non-adiabatic effects should cause os- 
cillations |y,y_7| , so it would not surprise if numer- 
ical work eventually confirms same (but only for 
a certain region of parameter space) . 

2. QUANTUM KINETIC EQUATIONS 

The dynamical variable of interest is the 1-body 
reduced density matrix p for the Va/^s subsystem. 
There is a corresponding density matrix for the 
antineutrinos. These functions depend on time 
t, or equivalently temperature T, and momentum 
y = p/T. We will employ the decomposition 

p^^{Po + P-a). (6) 

The diagonal entries are conveniently normalised 
distribution functions /„ and /s for Va and i^^, 
respectively: 

The off-diagonal entries Px.y are the coherences. 
The Quantum Kinetic Equations (QKEs) are 
II 

^2. (8) 



dP 



dPo 
dt 



P- 



)} 



(9) 



-^ = V^P-D{Px± + Pyy) 

The damping or decoherence function D is given 
by D = ^, with F being the total collision rate 
of the weak eigenstate neutrino of momentum y 
with the background plasma, T{y) ~ Ka G'^T^ y, 

where Kg — 1.27 and k^ = Kt — 0.92 (for rUe ~ 

T r^ TO^). 

The function f^q is the Fermi-Dirac distribu- 
tion, 

1 



feqiy) 



(10) 



1 + ev-t^" 

where jla = fJ-u^./T is a dimensionless chemical 
potential, and f°g is the /i^^ = case. 

The quantity V is given by 

V = /3x + Az, 

where 

sin 29o cos 20o 



/3 



4 



A 



Va, 



(11) 

(12) 



and £o = 2yT/Am'^ is the vacuum oscillation 
length. 

The V X P term describes coherent nonlinear 
MSW evolution, the D entry represents collisional 
decoherence, while R quantifies repopulation of 
the Va distribution from the background plasma. 

An overbar will designate corresponding an- 
tineutrino functions. The kinetic equations are 
identical except for the substitution L -^ —L. 

The QKEs together with a + s lepton number 
conservation supply 



dL„ 



dt 



2n 



P{Py - Py 



J Jo 



,y_J_eq_ 

27r2 



dy 



(13) 



as the asymmetry evolution equation. This is a 
redundant equation, but numerically very useful 
in removing the problem of taking the difference 
of two large numbers to obtain the asymmetry. It 
is nonlinear, because A and A depend on L^^ . 

3. ANALYTICAL AND PHYSICAL UN- 
DERSTANDING OF ASYMMETRY 
GROWTH 

The growth curves in Fig.n^ were obtained by 
numerically solving the QKEs in conjunction with 
Eq.(n3). In this section I review how the QKE 
"black box" may be opened and the basic fea- 
tures of these graphs analytically and physically 
understood. This will provide a partial reply to 
the claims of Ref.f|l9|. A more complete riposte 



can be found in Ref. |2 

3.1. Collision dominated epoch 

The collision rate F is proportional to T^, so 
at high enough temperatures collisions dominate 
the evolution. In this epoch it is legitimate to 
take R ~ 0, because fa does not depart from 
equilibrium form very strongly. See Ref. pl| for a 
discussion. The QKEs then reduce to 




In the instantaneous diagonal basis Q = UP this 
equation becomes 






(15) 



where /C^ is the diagonal matrix of eigenval- 
ues fci.2,3 of /C. In the adiabatic limit we set 
dlA~^ /dt = and Eg. ([l5|) may be formally solved 
to yield IHII 



P(i) ~ U-\t)e^o '^^^''u{0)P{0) 



(16) 



Under most circumstances, the eigenvalues 
take the form |l7| 



ki = kl 
with 
d = D 



lUJ, 



(3^D 



dP 



(17) 



2 '' 



\^ + fi"^ + ky,D + -B. (18) 



In the high T collision dominated regime 
and 



i:»2 + A2 + ^^ 



<i?, 



(19) 



(20) 



except at the very centre of the resonance. For 
full technical details concerning the resonance 
regime see Ref. [p^. When collisions dominate, 
exp(/ ki.2dt') -> 0, and Eq.(|l|) yields 



Pyit) 



a 



Pzit), 



so that Eg. ([l3|) becomes 



(21) 



(22) 



One also obtains a sterile neutrino production 
equation which I will not display. These equa- 
tions involve the distribution functions only - the 
coherences have been eliminated - so we have a 
BoUzmann limit. The root of this lies with the 
collisional damping by —d of the oscillatory be- 
haviour driven by w. 

Substituting for the various functions on the 
righthand side of Eq.(p3) one finally obtains, to 
leading order [p]J^,|2^ , 



dL^ 



2rp3 



^[x + ic-b + ay][x+ic~b-a)^y ^ ' 



where 



{DIq)"^. This has been called 



the "static approximation" . There are correction 
terms also which I will not discuss. If Arn^ > 0, 
then it is easy to see dL^^/dt ~ —L and so L^^ 
evolves so as to drive the effective asymmetry L 
to zero. The Am? < case is, however, quite 
different. The function b is then always positive, 
decreasing as T^, and the sign of the asymme- 
try derivative changes from negative to positive 
at a critical temperature Tc due to the evolving 
net effect of the (c — b) factor in the integrand. 
Immediately after Tc, a brief spurt of exponential 
growth occurs because now dL^^/dt '^ -\-L. 

The dashed line in Fig|2, taken from Ref.H, 
shows the asymmetry growth curve obtain from 
Eq. (p3|) . The dash-dotted curve displays asym- 
metry evolution as obtained from the full QKEs 
for the same oscillation parameter choice. The 
explosive growth at Tc is well approximated by 
Eq. (p3|) , demonstrating that collision dominated 
adiabatic evolution is the dominant physical ef- 
fect during this epoch. At lower T, the mag- 
nitude of the asymmetry is seriously underesti- 
mated. New physics has taken over B . 

3.2. MSW dominated epoch 

We begin by observing that Eq. (p3) reduces to 
dLu^/dt — when the collision rate F vanishes. 
This helps to explain why the approximate evo- 
lution equation reviewed above underestimates 
asymmetry growth at lower temperatures. When 
collisions are unimportant, the QKEs describe 
pure (non-linear) matter-affected evolution, so we 
now turn to a study of the MSW effect. 

Putting D = F/2 = in the adiabatic evolution 
matrix K. produces the eigenvalues 



ki=k*= iV/32 + A2, fca = 0. 



(24) 



dt 



Att^t 



There is no damping so the evolution is oscilla- 
tory, with |fci.2| being the matter-affected oscil- 
lation frequency for a neutrino of momentum y. 
The conditions for the validity of Eqs.(22) and 
(H) no longer hold 0,||. 

Furthermore, the condition of the plasma is 
now different because of the significant neutrino 
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Figure 2. Asymmetry growth curves from the QKEs (dash-dotted line), the static approxmation (dashed 
hue) and undamped adiabatic MSW transitions (sohd Hne) as presented in Ref.[Q] for Arri^ = — 10 eV^ 
and sin^ 26*0 = 10"^. 



asymmetry, leading to a separation of the neu- 
trino and antineutrino resonance momenta y^ [^ . 
Writing a = ^aU and b = £,by^, the MSW reso- 
nance conditions c — 6 =F a = imply that 



Vu 



and 



yh 



ia + 


Ve^ + 46 




26 


-Ca 


+ vea+^^b 



26 



(25) 



(26) 



If i > 0, then 6 is positive and y^ quickly moves 
to the high-momentum tail of the distribution. 
The partial cancellation between the terms in the 
numerator of y^ keeps the antineutrino resonance 
within the body of the distribution. Doing the 
algebra with 6 ^ 6 one obtains 



1 



y-'- 



8V2C{i)GFL.^T^' 



(27) 



where we have taken L ~ 2^^^. The resonance 
momentum evolution p] is shown in Fig.H. 



The separation of the resonances means that 
the MSW effect has very asymmetric (pardon the 
pun) consequences for neutrinos and antineutri- 
nos: neutrino oscillation are strongly suppressed, 
whereas 77q's get MSW converted into sterile 
states. This keeps the positive L growing Q. 
(The final sign of L, but not its magnitude, de- 
pends on the high T initial neutrino chemical po- 
tentials and so cannot be predicted.) The MSW 
effect is of course automatically taken care of 
by the QKEs. But in order to isolate its influ- 
ence, a more physical approach is useful. We 
will stay with the L > case. For small res- 
onance widths and 100% (adiabatic) conversion, 
the rate of change of the asymmetry is propor- 
tional to the speed with which the antineutrino 
resonance moves through the distribution 0: 



dL„ 



dT 



T3 



[hJy^)-hM)]iyv) 



2^ 

dT 



= -X 



dT 



(28) 
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Figure 3. The separation of the neutrino and antincutrino resonance momenta as the asymmetry grows. 
The oscillation parameters are as for Fig.|| This figure first appeared in Ref.[^. 



Using Eq. (E^) one then obtains 
dL,^ _ 4 XyL 



dT 



^1+S 



(29) 



as a non-linear evolution equation for the asym- 
metry, for the case where -rf < 0. We employ in- 
stantaneous repopulation before neutrino decou- 
pling, so fa is understood to be of Fermi-Dirac 
form with the appropriate (evolving) chemical po- 
tential. After chemical and kinetic decoupling one 
has to handle repopulation in a more complicated 
way: see Refs.[|[|2] for further explanations. 

The solid line in Fig|| shows the result of nu- 
merically integrating Eq. (p^) , beginning at T ~ 
Tc/2. The agreement with the QKE result is 
excellent, which shows that the MSW effect has 
taken over by T ~ Tc/2 Q. Furthermore, when 
L^^ <C 1, Eq. (p9|) reduces to 

^^ ~ -^^ (301 

which immediately implies that L^^ ^ T^'^. The 
power law behaviour displayed by the QKE solu- 



tion is thus easily understood as a consequence of 
undamped adiabatic MSW transitions p|. 

As the asymmetry continues to grow, the an- 
tincutrino resonance eventually moves out of the 
body into the high momentum end of the dis- 
tribution. The asymmetry becomes frozen at 
some final steady state value, whose approximate 
magnitude can be easily understood by integrat- 
ing the Fermi-Dirac distribution from y '--^ to 



L 



final 



1 



T„ 



4C(3) V T 



1 + ey 



3 /T 



8\n 



(31) 



The temperature ratio takes care of reheating due 
to e+e" annihilations at T ~ me — 0.5 MeV. 
Equation ( pi|) gives the approximate magnitude 
only, because the distribution changes with time 
as the asymmetry is created. Numerically, the 
final values found by incorporating proper ther- 
malisation effects are quite close to this estimate. 



4. CONCLUSIONS 

Active-sterile neutrino oscillations will produce 
a final steady state asymmetry of order 3/8 pro- 
vided the oscillation parameters are in the appro- 
priate range. Applications include the suppres- 
sion of sterile neutrino production prior to BBN 
and the alteration of the primoridal Helium abun- 
dance due to a z/g asymmetry. Inhomogeneous 
baryogenesis can seed lepton domain formation 
leading to inhomogeneous BBN. The asymmetry 
growth curves obtained from brute force numer- 
ical solution of the QKEs have been understood 
analytically and physically. 
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